Light-intensity ssusceptibility and "active" noise spectroscopy 
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In this paper, we consider informative potentialities of the "active" optical noise spectroscopy, 
under which we understand, generally, spectroscopy of response of a multilevel quantum system to 
the resonant optical field with its intensity modulated by "white" noise. We show that calculations 
of such a response can be most conveniently performed, in the linear approximation, by introducing 
the notion of light-intensity susceptibility (LIS) whose spectrum is determined by Laplace transform 
of the response to a small step-wise change of the optical field intensity. The results of calculations 
for a simple four-level quantum system show that its LIS spectrum may provide information not only 
about the ground-state structure (like conventional Faraday-rotation-based spin noise spectroscopy), 
but also about properties of the optical transitions (including nutation frequencies in the applied 
optical field). From the experimental point of view, such a noise spectroscopy of the intensity-related 
susceptibility can be especially efficient in combination with the up-to-date spectrum analyzers 
providing extremely fast data processing. 



I. INTRODUCTION 

Optical spectroscopy, with its frequencies lying in the 
range of 10 15 Hz, is known to be capable of studying spec- 
tral features of the system at frequencies many orders of 
magnitude lower (like, e.g., Zeeman or hyperfine split- 
tings). Such a high-resolution spectroscopy often em- 
ploys, for detection of these features, the intensity cor- 
relation characteristics of the light coupled to the sys- 
tem. Ideology of the light intensity noise spectroscopy 
has been primarily realized by Forrester [l| and Hunbury- 
Brown and Twiss @ and then has been further developed 
for studying dynamics of motion of different atomic and 
molecular systems including dynamics of nonstationary 
quantum states (see, e.g., [3|]). With the advent of lasers, 
the light intensity noise spectroscopy has crucially in- 
creased its sensitivity and gained a deeper practical sense. 

The methods of noise spectroscopy intended for study- 
ing spectra of eigen-frequencies of quantum systems can 
be divided into two main classes. In the measurements of 
the first class, optical probing of the sample is supposed 
to be nonperturbative: the light beam passing through 
the sample in the region of its transparency acquires addi- 
tional intensity or polarization fluctuations that carry the 
sought information. Ideally, intrinsic noise of the probe 
beam, in this approach, is negligibly small. A typical ex- 
ample of the experimental technique of this class is the 
Faraday-rotation-based spin noise spectroscopy (SNS), 
which provides a unique opportunity of studying mag- 
netic resonance and spin dynamics of atomic and solid- 
state systems practically in a perturbation- free way [IJ-0] . 
The recent enormous growth of interest to this experi- 
mental technique, primarily demonstrated more than 30 
years ago [1], is associated, to a considerable extent, with 
great progress in the present-day electronics, which has 
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made it possible to digitize arrays of data with the sam- 
pling rates up to several GHz and to perform fast Fourier 
transform (FFT) in real time. As a result, the new class 
of spectrum analyzers with FFT has allowed to increase 
sensitivity of the measurements and, correspondingly, to 
shorten the acquizition time by several orders of magni- 
tude (as compared to the standard sweeping spectral an- 
alyzers). Due to this technical refinement, SNS has been 
successfully applied to semiconductor systems, including 
microsamples and nanostructures [6] and is now about to 
turn into a routine tool of magnetic spectroscopy. 

The other experimental approach to the noise spec- 
troscopy that has been developed in 90s and applied to 
atomic systems, is based on the use of a randomly mod- 
ulated light beam acting on the system in the region of 
its absorption [8l-fl3T|. This technique is deliberately per- 
turbative: the system excited by the pump light modu- 
lated in a wide frequency range responds selectively to 
certain resonant frequencies, and, as a result, the noise 
spectrum of the transmitted light appears to be modi- 
fied. In the above papers, as a source of light was used 
a diode laser with specific correlation properties of its 
output emission: intensity fluctuations were extremely 
low (they were claimed to lie below the shot-noise level 
[l^|). whereas its frequency noise was well pronounced. 
The frequency modulated light, due to interaction with a 
narrow spectral line of the atomic medium, produced the 
required noisy excitation of the system. It was shown, 
in those experiments that the atomic system, under cer- 
tain conditions, may exhibit resonant response not only 
at the frequency of the ground-state magnetic resonance 
(Larmor precession), as it was primarily demonstrated 
by Bell and Blum|15|, but also at Rabi frequency of the 
resonant transition [10], which may substantially widen 
informative potential of the noise spectroscopy. 

It is noteworthy that this kind of noise spectroscopy 
was never considered as a possible experimental tool for 
studying a wider circle of materials, including semicon- 
ductor systems and nanostructures, highly important for 
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contemporary physics and applications. At the same 
time, from the methodological viewpoint, experimental 
measurements of the intensity-related susceptibility using 
fluctuating optical fields have much in common with the 
technique of conventional "passive" noise spectroscopy 
and thereby have all chances to take an advantage of 
the new instrumental opportunities provided by the up- 
to-date electronics. 

In this paper, we consider a general theoretical model 
that describes the effect of a four-level system on the in- 
tensity noise spectrum of the transmitted resonant light. 
The light intensity noise spectrum is assumed to be 
"white" and essentially exceeding in magnitude the shot- 
noise level. We analyze interaction of the light beam 
with the medium by introducing the notion of the light- 
intensity susceptibility (LIS), which in our opinion, may 
be useful for solving certain problems of nonlinear optics. 
In particular, this notion was, in fact, used in the analysis 
of dynamics of a saturable absorber in modulated optical 
fields HI]. 

In what follows we study spectral behavior of the 
intensity-related susceptibility in the range of relatively 
low frequencies (e.g., of the ESR range) . We will show 
that frequency dependence of the light-intensity suscep- 
tibility contains information not only about energy and 
relaxation characteristics of the system accessible for lin- 
ear spectroscopy, (including conventional "passive" spin 
noise spectroscopy), but also about specific dynamics of 
optical transitions (nutation frequencies) in the field of 
the noisy resonant pumping. We will see that, in contrast 
to the conventional linear optical spectroscopy, which al- 
lows one to measure frequencies of transitions between 
the levels of unperturbed states of the system, the LIS 
spectroscopy makes it possible (for sufficiently long re- 
laxation times) to measure frequencies of transitions be- 
tween the states of some effective Hamiltonian dependent 
on the optical field intensity io-03 



II. STARTING POINTS 

We define the light-intensity susceptibility (LIS) in the 
following way. Assume that the system under study 
is subjected to action of a quasi-monochromatic opti- 
cal field of the frequency oj, with its intensity / hav- 
ing a dc (Iq) and weak ac (8I(t)) components / = 
Jo + 8I(t). Correspondingly, the power P absorbed by 
the system will also have the dc and ac components: 
P = Po + SP(t). We define the light-intensity susceptibil- 
ity as a linear susceptibility x( l connecting the ac com- 
ponents of the intensity and absorption in the case when 
81{t) ->■ 6Ie- % &,8P(t) ->■ SPe-^ and SP = x«)SI- 

A possible experimental observation of frequency de- 
pendence (spectrum) of the light-intensity susceptibility 
with the use of the spectrum analyzer looks as follows 
[Fig. []Ja)]. The laser beam subjected to a broadband 
noise intensity modulation [18j passes through the sam- 
ple under study and hits the photodetector, whose out- 



M 



|t3/2> - 



CP 



D 

■Or 



|t1/2> 



NG 



SA 



Figure 1. (a) schematic of experimental setup for observation 
of light intensity susceptibility. L - laser, M - optical modu- 
lator, CP - circular polarizer, S - sample, D - photodetector, 
SA - spectrum analyzer, NG - noise generator, (b) energy 
structure of the model system. 



put signal is fed to the input of the spectrum analyzer. 
If the intensity modulation spectrum of the laser beam 
is "white" , then we will observe, at the output of the 
spectrum analyzer, the LIS spectrum of the studied sarn- 
ie. The experiments of this kind have been described in 
[T]| where the resonant paramagnetic susceptibility of 
atomic gases has been detected. It should be emphasized 
that, though we consider a linear intensity susceptibility, 
it can be revealed only by optically nonlinear media. In 
the next section, we will calculate the LIS of the simplest 
model system and will study its spectrum. 



III. MODEL CALCULATIONS 

The sample under study is assumed to be an ensemble 
of particles, and the total power absorbed by the sam- 
ple irradiated by the electromagnetic field to be equal 
to the sum of contributions of individual particles. The 
energy spectrum of individual particle is assumed to con- 
sist of two doublets [Fig. QJb)]: the ground | ± 1/2) and 
the excited | ± 3/2) one (here we indicate projection of 
the angular moment onto the light propagation direction, 
which is assumed to coincide with the z-axis of the our 
coordinate system). This energy structure arises, in par- 
ticular, in the atomic multiplet perturbed by an axially 
symmetric field. [19] In this case, magnetic splitting of 
the ground doublet is characterized by isotropic g-factor, 
whereas the doublet |±3/2) is split only by the field paral- 
lel to the z-axis. Spectral distance between the doublets 
(denoted by f2) is supposed to be close to frequency to of 
the optical field cj rj f2 and essentially higher than other 
frequency parameters of the problem. 

We will calculate the LIS spectrum we are interest- 
ing in as follows. Let us assume that the system under 
study was subjected, during sufficiently long time, to a 
monochromatic optical field with fixed intensity, so that 
the power absorbed by the system has already reached 
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its steady-state value. Let the amplitude of this field now 
experience a jump that induces a transient process cor- 
responding to transition to a new stationary state. Be- 
low, we will find the transient characteristics (i.e., tempo- 
ral behavior of the absorbed power during this process), 
which is connected with the sought susceptibility by the 
Laplace transform. 

Let us pass to implementation of the above program for 
the conditions of the experiment [la ]. The optical field 
acting upon the system is assumed to be circularly po- 
larized and the whole system to be placed in a transverse 
magnetic field H directed along the x-axis [Fig.QJa)]. 

Since the magnetic splitting of the ground doublet 
is supposed to be small compared with the optical fre- 
quency u>, among the matrix elements of the operator of 
interaction of the system with optical field (denote it Hi), 
we have to take into account only those corresponding to 
quasi-resonant transitions: 



We will mark with tilde all the operators having the 
form H = SHS. Then, direct calculations show that 
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And the density matrix p = SpS meets the equation 
idp/dt = [H, p]. Transition to the "rotating coordinate 
frame" is performed using the following time-dependent 
transformation of the operators. Let us pass to the hatted 
operators and new density matrix using the rule 



(-l/2|ffi| 
(+l/2\Hi\ 
(+3/2|#i| 
(-3/2|ffiJ 
<+l/2|ffi| 
<-l/2|#i| 



+ 1/2) 
-1/2) 

+ 1/2) 
-1/2) 
+ 3/2) 
-3/2) 



non resonant 
non resonant 

Ae~* ut 
Ae lujt 
Ae lujt 
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Here, the constant A describes intensity of the transitions 
under circularly polarized excitation. If we arrange pos- 
sible states of the system in the order | + 1/2), | — 1/2), | — 
3/2), | + 3/2), then the matrix Hi will have the form 



Hi = A 



/0 fj\ 

0/20 

p, 

\fi 0/ 



fx = e lut 



(1) 



At A = 0, the Hamiltonian matrix of the system, in the 
basis introduced above, will have the form 



'0 v N 

v 

ft 

vO ft/ 



2h 



(2) 



At A ^ 0, the Hamiltonian of the system depends on 
time and equals H = Hq + H\. Neglecting the relaxation 
processes (they will be taken into account later) , dynam- 
ics of the density matrix of the system p can be described 
by the Liouville equation ip = [H, p\. To treat this equa- 
tion, it is convenient to pass to the basis, in which Hq 
is diagonal, and then to the "rotating coordinate frame", 
where the total Hamiltonian is time-independent. The 
first step is made using a unitary transformation with 
the matrix 



5 = -r 



(11 \ 

1-10 

y/2 

\0 V2~J 



with S" 1 = S (3) 



H = e lMt He- lMt o = e lMt pe- lMt (6) 
Where the matrix M has the following elements 



M = 



/0 0^ 



oj 

\0 ojj 



(7) 



Then, one can easily see that a satisfies the equation 
la = [H — M, a) . Calculations of the matrix H show 
that it contains time-independent elements and elements 
proportional to e ±2lujt . The matrix elements oscillating 
at double frequency are ignored as essentially nonreso- 
nant and not affecting dynamics of the density matrix 
a. Thus, the matrix H, with allowance for the aforesaid, 
can be written in the form: 



H 



>v ^ 

-v d 

ft 

yd d ft/ 



A_ 
V2 



(8) 



It follows from this expression that the third column and 
third row of the matrix H are mutually diagonal and can 
be omitted. Now, all the operators can be presented by 
3x3 matrices, and the equation of motion for the matrix 
a acquires the form 







i& = [W,a], 



W 



A = ft-i 



(9) 

with the time-independent matrix W. 

As mentioned above, our task is to calculate transient 
dynamics of the power P absorbed by the system after a 
jump of the optical field amplitude. This power can be 
expressed in the form 



P 



dt 



Sp pH 



dt 



Sp aH = Sp &H- 



Sp.f (10) 
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The last term in (fTU)) can be neglected because the op- 
erator H oscillates in time at double optical frequency, 
and the mean value of the corresponding contribution 
vanishes. Thus, 

P = Sp u PL = — i Sp [H - M, a]H = i Sp [M, a]H 

(11) 

Using explicit formulas (8) and (9) for the matrices M 
and H, we finally have 



P = 2ujd 



Im CT13 + Im (T 2 3 



(12) 



Note that account of relaxation terms in the equation 
of motion for the density matrix leads to appearance of 
terms proportional to the relaxation rates in the expres- 
sion for the absorbed power P. We will assume that 
these rates are much lower than the optical frequency oj 
entering Eq (fl~2|) and will use Eq. (fl~2|) in the presence of 
relaxation processes also. 

Now, we will write down Eq. (J5|) in components of the 
density matrix and introduce relaxation needed to obtain 
the steady-state regime of absorption: 



ion 
i&i 3 

l&U 
^33 



iT, 



CT12 [2v - iT± 2 ] + d[a 32 - CT13] 
-(72i [2v + iT^\ + d[a 3 i - CT23] 
(7i3 [v - A - iT^] + d[a 33 - (7n 
a 3 i [A - v + 
(7 23 [- A - v 

(7 32 [A + V + 

4°3i - 0-13] 

d[a 32 — (723] 
d[a 13 - CT31 - 



C12J 

31 j T d[crxi + (7 2 i - C33] 

lT 23~] + 4°33 - 0-21 - CT 22 ] 

CT ii) T r 



32 j T d[(7 22 - 
1 +(ff22 



r (733 Ti 

H 033T]" 

023 — cr 32 



1 - ((722 - CTlljTf 1 



(is) 

The form of the relaxation terms for diagonal elements of 
the density matrix implies that states 1 and 2 are the low- 
est, and the temperature is so high that they are equally 
populated. State 3, on the contrary, is assumed to lie 



so high that its equilibrium population is zero. These 
assumptions correspond to the conditions g(3H << kT 
and Ml » kT. Let us introduce the following notations 
for the real and imaginary parts of the matrix elements 
a and relaxation times 
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(14b) 



(14c) 



(14d) 



Here, we denoted the dephasing time in the magnetic 
doublet as T2 and the dephasing time in the transitions 
2-3 and 1-3, saturated by the optical field, as T2. By 
defining the vector-column v as 



and the matrix G(d) as 
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Using Eq. (|13p . we obtain the following equation for v Since the last row of matrix G is the sum of two previ- 

ous ones, this matrix should always have a zero eigen- 

dv number. The corresponding eigen-vector (denote it e 9 ) 

Tt =G{d)v (17) 
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allows us to obtain the steady-state regime of the system 
saturated by the optical field. The vector b, correspond- 
ing to the steady-state matrix of the system, evidently 
represents the vector e 9 normalized as b = q~ 1 e 9 , where 
q = ej + e| + eg. In accordance with the formulation 
of the problem, given above, we assume that the density 
matrix of the system, before the jump of the optical field 
amplitude, was determined by the vector b. In this case, 
the power P Q absorbed by the system before the jump 
can be expressed through the components of this vector 
using Eq. (TT2]) 



P = ud[b 5 + b 6 



(18) 



If the jump of the optical field amplitude d — > d + Sd 
occurs at t = 0, then, to find dynamics of the system after 
the jump, we have to solve the equation v = G(d + Sd)v 
with the initial condition v(t = 0) = b. Denote by p l , i — 
1,2,..., 9 the eigen-vectors of the matrix G(d + 5d) and 
decompose the vector b over them: 



E^ 1 



(19) 



Vectors p l and coefficients Ci can be easily found nu- 
merically. Now, the vector v describing dynamics of the 
system after the jump can be written in the form 



iP e 



(20) 



Where Ai are the eigen-numbers of the matrix G(d + Sd). 
Dynamics of the absorbed power is given by the equation 



P(t) = u(d+8d)yC i e* t lp\ 
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Pel 



(21) 



The quantity SP = P(t) — Pq can be considered as 
a response of the system to the jump 5d, and, at small 
5d, this response should be linear. The LIS we are in- 
terested in, is evidently proportional to the susceptibility 
K(i^) connecting SP and Sd. Then, if S(£) is the Laplace 
transform of the response to the jump, 



s(0 



lim 
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p oo " 

/ e lit [P(t)-P ]e- at dt= -oj{d + Sd) V 
Jo „■_-, 



it, + Xi 



ujd(b 5 + b 6 ) 



(22) 



then \K{i£)\ 2 — \£,S{£)\ 2 , and, for the LIS spectrum, we 
obtain 

lx(*0l 2 ~l£S(0l 2 ( 23 ) 

Thus, calculations of the LIS spectrum, within the sim- 
plest model described above can be performed using Eqs. 
(gH and d23J|. 

IV. RESULTS 

Let us calculate the LIS spectrum for essentially dif- 
ferent relaxation rates in the optical transitions. 

When the relaxation times entering the equations for 
elements of the density matrix (TTH1) substantially exceed 
the period of optical nutation d, the dynamics of the den- 
sity matrix for times smaller than the relaxation times is 
determined by Eq. ©. Under these conditions, this dy- 
namics represents oscillations at frequencies equal to all 
possible differences between the eigen-values of the effec- 
tive Hamiltonian W. The quantities Wi,i — 1,2,3, in 
this case, can be calculated either numerically or using 
Cardano's formula. Thus, the LIS spectrum will reveal 
peaks at appropriate frequencies, which, in this case, will 
be three: £ x = |Wi-W 2 |,6 = |Wi-W s |,6 = \W 2 -W 3 \. 
This is demonstrated in Fig. [Ha), which shows the LIS 
spectrum obtained using Eqs. (|2"2")l and ([2"B")) for d^ 1 << 
T\,Ti,T2. Verical lines indicate frequencies — 1, 2, 3. 



In the absence of the magnetic field, one can easily see 
that W h2 = [A ± VA 2 +4d 2 ]/2, ' W 3 = 0. At 
d/A << 1, the LIS spectrum shows peaks at frequencies 
£ w A and £ w d 2 / A (the peak is weakly pronounced), 
while at d/A << 1, at the frequency of optical nutation 
2d and at d (the peak is also weak). 

The situation strongly changes when passing to short 
relaxation times of the optical transitions T2, ti << d . 
This case corresponds to the experiment [15] when the 
system becomes relatively transparent for the circularly 
polarized light with its intensity oscillating at the fre- 
quency of magnetic splitting (ESR frequency). In this 
case, the LIS spectrum shows a dip at the frequency 
of magnetic splitting 2v [Fig. [2Kb)] . The condition 
72, ti << d^ 1 allows us to neglect the time derivatives 
<7i3, 1731, o"23, 032, 033 in f| 1 3 [) as compared with the relax- 
ation terms. This, in turn, makes it possible to express 
algebraically the matrix elements 013, 031, 023, 032, C33, 
pertinent to optical transitions through the elements 
C12, C21, fTn, CT22, related to the magnetically split dou- 
blet and, thus, to obtain the relationships similar to those 
presented in [15j . 



V. CONCLUSION 

The goal of this paper is to show, using as an example 
a four-level system, potential usefulness and efficiency of 
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Figure 2. (a) the LIS spectrum for relatively long relaxation 
times in the optical transitions T2 = 10, ri = 2t2, (b) the 
same for relatively short relaxation times T2 = 0.01, ri = 2t2. 
The rest parameters are the same for both figures and are: 
T 2 = 10, Ti =2T 2 ,d = l,v = A = 1.5. 

the "active" noise spectroscopy that employs transfor- 
mation of intensity noise spectrum of the light passing 
through the medium as a source of information about its 
eigen-frequencies. This type of the noise spectroscopy 
has much in common with the conventional Faraday- 
rotation-based spin noise spectroscopy and, under cer- 
tain conditions, provides the same information about the 
system and also may serve as a method of magnetic spec- 
troscopy. At the same time, the technique considered 



here, though deprived of the advantage of being non- 
perturbative, has a number of specific properties and un- 
questionable merits. In contrast to the Faraday-rotation- 
based technique, this method allows one to get informa- 
tion about optical transitions and properties of the ex- 
cited states, rather than only about the ground-state dy- 
namics and structure. The signal of the "active" noise 
spectroscopy, being a result of induced (rather than spon- 
taneous) response of the system, can be controlled by the 
intensity of the optical field and, evidently, can be easier 
detected. This spectroscopy, unlike the "passive" FR- 
based, implies pure intensity-related measurements and, 
in the proposed version, does not require any polarimet- 
ric sensitivity. This technique employs fluctuations of 
the optical field rather than thermodynamic fluctuations 
of the ensemble of particles, and, therefore, the detected 
signal is not connected explicitely with the size of the 
light spot. It may seem strange, but a serious problem 
for this kind of spectroscopy may be the problem of laser 
source with a sufficiently high level of intensity noise in 
a sufficiently broad frequency range. Perhaps, special ef- 
forts will be needed to create a special laser source for 
the spectrometer of this kind. In principle, such a spec- 
troscopy of the light intensity susceptibility can be imple- 
mented with harmonically (rather than randomly) mod- 
ulated light. But in this case, to obtain the spectrum, 
we will have to scan the modulation frequency, and the 
experimental technique will lose its main advantage. 

We consider the spectroscopy of the light intensity 
susceptibility to be a promising experimental technique 
for studying energy structure and dynamics of multilevel 
quantum systems under resonant optical excitation. 
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